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ABSTRACT 

As the 50th anniversary of the Cerro Tololo Inter- American Observatory (CTIO) draws near, the 
author was surprised to learn that the published latitude and longitude for CTIO in the Astronomical 
Almanac and iraf observatory database appears to differ from modern GPS-measured geodetic posi- 
tions by nearly a kilometer. Surely, the position for CTIO could not be in error after five decades? 
The source of the discrepancy appears to be due to the ~30" difference between the astronomical and 
geodetic positions - a systematic effect due to vertical deflection first reported by Harrington, Mintz 
Blanco, & Blanco (1972). Since the astronomical position is not necessarily the desired quantity for 
some calculations, and since the number of facilities on Cerro Tololo and neighboring Cerro Pachon 
has grown considerably over the years, I decided to measure accurate geodetic positions for all of the 
observatories and some select landmarks on the two peaks using GPS and Google Earth. Both sets 
of measurements were inter-compared, and externally compared to a high accuracy geodetic position 
for a NASA Space Geodesy Program survey monument on Tololo. I conclude that Google Earth can 
currently be used to determine absolute geodetic positions (i.e. compared to GPS) accurate to roughly 
±0.15" (±5 m) in latitude and longitude without correction, or approximately ±0".10 (±3 m) with 
correction. I tabulate final geodetic and geocentric positions on the WGS 84 coordinate system for 
all astronomical observatories on Cerro Tololo and Cerro Pachon with accuracy ±0".l (±3 m). One 
surprise is that an oft-cited position for LSST is in error by 9.4 km and the quoted elevation is in 
error by 500 m. 

Subject headings: astronomical databases: miscellaneous - reference systems - telescopes - reference 
systems 



1. MOTIVATION 

For 50 years, the Association of Universities for Re- 
search in Astronomy (AURA) 3 has supported a growing 
numbe r of ast ronomical observatories on ICerro Tololo 
and Cerro Pachon. These observatories include ones op- 
erated by the National Optical Astronomy Observatories 
(NOAO) 4 and several other "tenant" scientific research 
facilities. Since the founding of the Cerro Tololo Inter- 
American Observatory and purchase of the "El Totoral" 
property 25 November 1962 5 , CTIO has provided a crit- 
ical platform for astronomical observations of the south- 
ern skies 6 , and continues to support impressive science 
on no less than two dozen facilities. The neighboring 
peak Cerro Pachon has been developed over the past two 



Electronic address: emamajek@pas.rochcster.edu 

1 Cerro Tololo Inter- American Observatory, Casilla 603, La Ser- 
ena, Chile and National Optical Astronomy Observatory, 950 N. 
Cherry Ave., Tucson, AZ, 85721 

2 Current address: University of Rochester, Department of 
Physi cs fc Astronomy, Rochester, NY 14627-0171, USA 

3 I http : / /www .aura- astronomy, org / 1 

4 http://www.ctio.noao.edu/noao/ and http://www.noao.edu/ 

5 "Brief History ot the (Jerro Tololo Inter- American Ub- 
servatory" writt en by Victor Blanco in February 1993: 
|http: / /www. ctio.noao.edu/noao/content / ctio-history| 

u Including, of course, support for observations with the Blanco 
4-m telescope for the High-Z Supernova Search team and Super- 
nova Cosmol ogy Project which led to the 2011 Nobel Prize in 
physics (e.g . |Riess et aT]|1998| |Schmidt et al.||1998| |Perlmutter| 
|et al.|[T999| , and i mportant pioneering work by the Calan/Tololo 
KlNe la survey (e.g. |Hamuy et aJ.|1996p ■ CTIO continues to pro- 
vide support for the investigation ot tne cosmological acceleration 
through the Dark Energy Survey (DES) and Dark Energy Camera 
(DECam) project on the Blanco 4-m telescope. 



decades, hosting the SOAR 7 and Gemini South 8 tele- 
scopes. The Large Synoptic Survey Telescope (LSST) 9 
is currently under construction on the El Pen on peak of] 
Cerro Pachon. With excellent seeing, weather, monitor- 
ing and mitigation of light pollution, and strong legacy 
of development of infrastructure and trained staff, the 
Cerro Tololo and Cerro Pachon sites are well positioned 
to host astronomical observatories for the next half cen- 
tury. 

The precise positions of the observatories on Cerro 
Tololo has been historically useful for a variety of pur- 
poses, including lunar occ ultations (e.g. Lasker et al 



|1973| |Vilas fc L askcr 1 9771), planetary o ccultations [e.g. 
|Frencn~et al.||l 983; Hubbar d et al.|[l997 ), planetary ring 



occultati ons (e.g. [E lliot et al .|l981 ), and asteroid astrom- 
etry (e.g. Buie et al.||2U12| ). Among other purposes, pre- 
cise positions for telescopes on these peaks are important 
for applying velocity corrections for precise radial veloc- 
ity measurements with the CHIRON spectrogra ph on the 
SMARTS 1.5-m telescope ( |Schwab et aTTpHOl ), and will 
be important for accounting tor geocentric parallax when 
determining orbits of nearby small solar system bodies 
imaged with the Dark Energy Camera (DECam) on the 
Blanco 4-m telescope and LSST. 

Given the importance of a precise position for astro- 
nomical calculations, it came as a surprise to the au- 
thor that the position for CTIO published in the As- 
tronomical Almanac and in iraf obsdb.dat file (obser- 

http://www.soartelcscope.org/ 
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vatory database) differed from GPS-measured geodetic 
positions on Cerro Tololo by approximately a kilometer. 
As I discuss later, this is not an error - the discrepancy 
appears to be due to the mismatch between the astro- 



nomical and geodetic coordinates for Cer ro Tololo ( Har- 
rington et al.|1972}|Blanco fc Mayall|1972 ). However, the 



issue raised awareness that (1) there are conflicting (and 
sometimes incorrect) coordinates and elevations for fa- 
cilities on Cerro Tololo and Cerro Paction on the WWW 
and in the literature, and (2) coordinates have not yet 
been estimated for many of the new facilities on these 
peaks. 

This document is a summary of my notes and obser- 
vations on the subject of the positions of facilities on 
Cerro Tololo and Cerro Pachon. It may be updated later 
if improved measurements or corrections come to light. 
After some historical review and document archaeology, 
I present new measurements of accurate geodetic and 
geocentric coordinates for the astronomical facilities on 
Cerro Tololo and Cerro Pachon. Readers may prefer to 
skip ahead to Table [6] to see the final compilation of po- 
sitions. 

2. BACKGROUND ON COORDINATE AND ELEVATION 
SYSTEMS 

Following the ISO 6709 standard, geodetic positions 
in this paper are quoted as latitude (0), longitude (A), 
and elevation (H) on the World Geodetic System 84 
(WGS 84) coordinate frame, unless otherwise noted. El- 
evations will be discussed later in this section, as they 
require a more detailed explanation. North latitude and 
east longitude are positive. Longitude 0° on the WGS 84 
frame corresponds to the International Reference Merid- 
ian defined by the International Earth Rotation and Ref- 
erence Systems Service (IERS), which lies approximately 
5" (~100 m) east of the meridian at the Royal Observa- 
tory in Greenwich. One degree (°) of latitude <fi is 110.86 

km, one minute ( ) is 1.85 km, and one second ( ) is 30.8 
m. At the latitude of Cerro Tololo, one degree of longi- 
tude A is 96.32 km, one minute ( ) is 1.61 km, one second 
is 26.8 m. 

Global Position System (GPS; also called NAVSTAR) 
is the well known navigation system supported by a con- 
stellation of 24 orbiting satel lites launched by the U.S 
Depa rtment of Defense (e.g. Hofmann-Wellenhof et al 
1992|. Altitudes measured by GPS are referenced to a 



gravitational equipotential surface that defines sea level: 
the World Geodetic System 84 (WGS 84). According 
to a Department of Defense publication tilted "Mili- 
tary Standard for Department of Defense World Geode- 
tic System (WGS)" 10 , the WGS 84 provides "the ba- 
sic reference frame (coordinate system), geometric fig- 
ure for the earth (ellipsoid), earth gravitational model, 
and means to relate positions on various geodetic datums 
and systems for DoD operations and applications. 11 The 
WGS 84 geometric figure is an oblate spheroid with ra- 
dius at the equator of 6378137 m, and a flattening of / 
= 1/298.257223563 (with a corresponding polar radius of 
approximately 6356752 m. While there is a 1984 version 
of the geoid model, there have been later updates of the 



11 



Document MIL-STD-2401, dated 

J anuary. 1994, |http://earth-| 



Earth Gravitational Model (EGM; e.g. EGM96) - in- 
deed improving the resolution and accuracy of the geoid 
model is a never-ending pursuit of the geodesy commu- 
nity 12 . 

Quoted latitudes are usually geodetic ((/>), which mea- 
sure the angle between a normal at a position on a 
spheroid, and the equatorial plane. Sometimes geocentric 
latitudes (<f> ) are quoted, measuring the angle between 
the line between the surface position and the Earth's cen- 
ter, and the equatorial plane. Geodetic and geocentric 
longitudes A are identical as they share the same axis and 
reference meridian. Geodetic latitudes ((/>) and geocen- 
tric latitudes (4> ) at Tololo and Pachon differ by ~10'. 
Conversions for geodetic and geocentric coordinates can 
be found in Sec. 4 of the Explanatory Sup plement to the 
Astronomical Almanac ( Seidelmann 1992 1 . 



One can define geodetic height with respect to the ref- 
erence ellipsoid. Elevations are classically defined with 
respect to mean sea level, however this does not follow 
the ellipsoid exactly due to the local concentration of 
mass - and of course defining sea level on land is not 
trivial. Earth gravitational models (EGMs) are mathe- 
matical approximations for the geoid, the Earth's grav- 
itational equipotential surface. The gravitational vector 
is perpendicular to the geoid. The ocean's mean sea level 
roughly follows the geoid, so one needs to define the geoid 
in order to quote altitudes with respect to "mean sea 
level" on land, as is commonly done. The geoid varies 
from the ellipsoid shape by deviations of up to approx- 
imately ±100 m. In reality, the geoid is a very compli- 
cated and irregular surface following local concentrations 
of mass. In practice, it is defined by high order spheri- 
cal harmonic expansion expressions. One commonly used 
EGM is Earth Gravitational Model 1996 (EGM96). This 
defines the geoid with a spherical harmonic series of n 
— 360 (~100 km resolution). There are three different 
"heights" of note: geodetic height (h), orthometric height 
or elevation (H), and geoid height (N). They are related 



as: 



H = h — N 



(1) 



i.e. the elevation H is equal to the geodetic height h 
minus the geoid height N (measure d with respect to t he 
reference ellipsoid; see Chapter 4 of Sei delmarm||1992 ). 

An online tool from UNAVCO 13 was used to estimate 
the geoid height N on the EGM96 model for Cerro Tololo 
and Cerro Pachon. The mean geoid heights N for Tololo 
and Pachon are approximately 34.6 and 35.0 meters, re- 
spectively. Geoid heights at individual observatories will 
be tabulated at the end of the paper. EGM96 geoid 
heights across both peaks vary at the <0.1 m level from 
site to site. There are subtle differences in geoid models 
as they have improved in accuracy and resolution over 
the years. For example, at the position of the Blanco 4- 



11 For further details on WGS 
ations over the past decades 



and its minor alter- 
see NIMA Technical Report 



TR8350.2, "Department of Defense World Geodetic System 
1984, Its Definition and Relationships With Local Geodetic Sys- 
tems", 3rd edition, Amendment 1 (3 Jan 2000): http://earth- 
info.nsa.mil/GandG/publications/tr8350.2/tr8350_2. html 



info.nga.mil/publications / specs / printed/2401 / 2401 .pdt 



±-a http://www.iag-aig.org/ 

1 " http://www.unavco.org/community_science/science- 



support/geoid/gcoid.html 
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m telescope, the (EGM84, EGM96, and EGM2008 14 15 ) 
Earth gravitational models predict geoid heights N — 
(31.41, 34.62, 33.31) meters. Hence, the definition of 
"mean sea level" under the observatories has varied at 
the ±few meter level over the past decades. Through- 
out this paper, all geoid heights N refer to the EGM96 
model, unless otherwise noted. 

To what degree are Cerro Tololo and Cerro Pachon 
moving due to plate tectonics? While the common per- 
ception is that Chile is moving westward with the South 
American tectonic plate, there are surprisingly large vari- 
ations in the plate motion in Chile. These large varia- 
tions are traced through time series GPS measurements 
of ground stations in Chile maintained by the Jet Propul- 
sion Laboratory 16 . Generally, stations north of Santi- 
ago (Santiago, Valparaiso, Copiapo, Iqueque) are moving 
northeast with respect to WGS 84 at ~30 mm/yr, while 
Chilean stations south of Santiago are moving northwest 
at ^90 mm/yr. Although I have been unable to find 
accurate geodetic motions for Cerro Tololo and Cerro 
Pachon in particular, judging by the measurements for 
other well-studied benchmarks in Chile, it is very likely 
that the peaks are moving at <100 mm/yr, and most 
likely <30 mm/yr. UNAVCO 17 provides a plate mo- 
tion calculator (model GSRM vl.2) on their website, and 
at Cerro Tololo's position, they predict that the South 
American plate is moving with respect to the WGS 84 
coordinate system at 9.24 mm/yr (9.17 mm/yr N, 1.10 
mm/yr E) - however this velocity is substantially smaller 
than those reported for northern Chile JPL stations. 
Hence, over the five decade history of the observatory, 
the site has very likely moved <5 m (and indeed the UN- 
AVCO plate motion calculator would predict only ^0.5 
m of motion over 50 years). In summary, any deviations 
in the published positions for the observatory larger than 
these amounts are unlikely to be attributable to tectonic 
motion, especially on short timescales. 



3. PREVIOUS COORDINATES AND ELEVATIONS 

3.1. Coordinates 

Past published positions for Cerro Tololo and Cerro 
Pachon and individual structures are listed in Table Q] 
The Astronomical Almanac (A A; e.g. |U. S. Government 
Printing Office (USGPO)|2013| ) often listed positions for 
individual telescopes at major observatories during the 
early 1980 's, but only listed positions for 4 telescopes 
on Tololo in the 1982 edition: the Blanco 4-m reflector, 
the 1.5-m reflector, the 1-m reflector, and the 24" Cur- 
tis Schmidt telescope. These positions were pres umably 
adopted from th e CTIO facilities manual (e.g. I Walker | 



& Mufioz 19801). The AA ceased listing positions for 
individual telescopes at observatories after the mid-80 's, 
and provided only mean positions for observatories there- 
after. Note that the AA concedes that the observatory 
positions are a mix of geodetic and astronomical values, 
and the type of position is not provided. 

14 http:/ /earth-info. nga.mil/GandG/wgs84/gravitymod/egm20C 8/ 
|indexTritml| 



The Minor Planet Center (MPC) lists 18 two observa- 
tory code entries for CTIO: #807 (Cerro Tololo Obser- 
vatory, La Serena) and #102 (Cerro Tololo, La Serena- 
2MASS). At the time of writing (16 July 2012), the MPC 
lists identical longitudes A, and parallax constants pcos(f> 
and psiiuj) for these two codes: A = 289°. 1941 (degrees 
east of Greenwich), and pcos0 = 0.86560 and pshxj) = 
-0.49980. u (j>" in this case is apparently geocentric lat- 
itude, and should be labeled <f> to distinguish it from 
geodetic latitude <p. From these last two quantities we 
can calculate the parallax parameter p to be 0.999532, 
and latitude <j) = -30°. 00228 = -30°00'08". This is in 
good agreement with values derived later in this paper. 

3.2. 1972 Harrington et al. Determination of 
Astronomical and Geodetic Coordinates 



In an annual report, Blanco & Mayall (19721 reported 
the following regarding the coordinates for CTIO: 

11 A program to determine the precise geodetic and geo- 
graphic positions of Cerro Tololo was initiated by Dr. R. 
Harrington, U.S. Naval Observatory, Washington D.C., 
and Dr. and Mrs. V. M. Blanco, CTIO. Observations 
of the CEOS B satellite, which was flashed especially 
for this purpose, were made with the Curtis Schmidt 
telescope. The final reductions made at the Goddard 
Space Flight Center yield the following geodetic position: 
|0 = -3Cr > 10' 8". 2, \ = WW ; 48' 21" J] The geographic 
or astronomical position was determined by the method 
of equal altitudes with data obtained from a series of 
theodolite observati ons. The pr eliminary results are p] 
\= -SOP 9' 53" ±1" .8 A = WW 48' 54"±3"\ These 
figures suggest a deflection of the vertical of approxi- 
mately 39 arc sec in the west-northwesterly direction, 
approximately perpendicular to the orientation of the 
Andean Cordillera and the deep off-shore Chile-Peru 
oceanic trench, which are the probable sources of the 
deflection. " 



Examination of the original |Harrington et al. ( 1972 1 
study 19 suggests a typo in the deflection of vertical 
measured. The astronomical and geodetic coordinates 
for the Curtis Schmidt telescope on Tololo are listed in 
Table [l] |Harrington et al.| ( |1972[ ) concluded: 



http://gcographiclib.sourccforge.net/cgi-bin/GeoidEval 

16 http: / /sideshow. jpl.nasa.gov/mbh/scrics.html 

17 1 



http: //www. unavco.org/community .science/science- 



support / crustal_motion / dxdt / model.html 



"The astronomic coordinates of the Curtis Schmidt 
telescope can be compared to the above geodetic coor- 
dinates to obtain the deflection of the vertical. This 
deflection amounts to 30". H 1".7 towards an azimuth 
of 295° ± 3° , which corresponds to approximately 930 
meters on the ellipsoid surface. " 

Hence, the source of the discrepancy between the Al- 
manac and iraf positions for CTIO and what one would 
measure with a GPS or on Google Earth appears to be 
due to the difference in the type of latitude and longitude 
being reported (i.e. astronomical vs. geodetic). While 
the astronomical position is useful for calculating transit 

18 http: //www. minor planetcenter.net/iau / lists/O bsCodesF. html 

19 A copy of the Harrington et al. (1972) study proved 
exceedingly difficult to find, but the NUAU North library 

has a copy. A scanned copy can be downloaded at 

http: / /www. pas. rochester.edu/ ~emamajek/Harrington72.pdf 
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TABLE 1 

Published Coordinates for Cerro Tololo and Cerro Paction 



Site 



Ref. 



Elev. 
(m) 



Notes 



Cerro Tololo 



CTIO 

CTIO 



UT1U 



CTIO 



CTIO 



CTIO 



Mayall 



' (196 8) 
& Mayall 



Blanco' jTTCTayalI| (l972|| 
Harring ton ct ai. [| lb>Y2|l 
Walker k, Munoz [ pggu[ l 



Oliu 
CTIO 
Blanc o 4-m 
Bianc o 4-m 
Blanco 4-m 



SMAKlS 1 5-m 



SMAIti'S l.ti-m 



SMAitlS 1.3-m 
SMARTS 1.0-m 
SmaRTS u.y-m 
Schmid t u.b-m 
Schmid t U.b -m 
rTololo " O"' 

PltUMF 

LCOGT 



aa iyvb-iysu 

iraf 

AA 1983-2013 

MPC(2012) 

IWalker & Munoz 1 119801) 
AA 1982 . ' — 
iHubbard et al.l 119971 

AA 1982 ' *— — 1 

IWalker & Munoz 1 119801) 
MPCC2U12 1 — 



HilUlU \_/ j 

vlASA M onumen t 7401 1 

TiOIvIPT] 



Walker fc Munoz] (l980| l 
Walker fc Munoz ||l980| 
Harrin gton ct ai. " IlijY^ ll 
Harrington et aU (liyy^ll 
AURA survey notes 
NASA 

PROMPT website 
LCOGT website 
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AA 2013 
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Gemini website 
Zombcck (2007) 



SOAR 



SOAK 
SOAR 
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IjSST wepsite" 
LSST website 
Simms et al.| (|2005| ) 
Zombeckj ||2uU7| 
SOAR website 



-30 13 42 
-30 13 42 
-30 14 26.700 
-30 14 ... 
-30 10 20.1 



-30 14 16.8 
-30 21 ... 



-70 43 24 
-70 43 24 
-70 44 12.096 
-70 43 ... 
-70 48 00.1 



-70 44 01.4 
-70 49 ... 



2725 
2737 
2722 
2715 
2123 
2647 
2647 
2738 
2701 
2701 
2724.60 



Pachon monument (near SOAR) AURA survey notes -30 14 13.393 -70 44 04.396 

Note. — I list longitude as m easured with east being positive, following the 2013 Astronomical Almanac ^U. S. Gov-| 
|ernment Pri nting Offi ce (USGPO)| |2013 | ■ In the first column, hyperlinks are provided to plot the coor dinates on Google 
Maps^ for comparison. Notes: (q) |Mayall| jl968[ position is quoted longitude as A h 43 m 16 s W. (b) | B lanco ^^^^^^aH 



^^72 l geode tic position determined trom observations of GEOS B satellite with Curtis Schmidt telescope trom~Paj^^^m 
et al.| <1972| stu dy. |Blanco fc Mayall| l |1972 I also lists the [Harrington et al"| ( | 1 9 astronomical coordinates. (c]~^^rrin^ 
al.l fll972|| "approximate coordinates, taken from a topographic sheet 6i 



ton et al.| | [ 

lipsoid and tti 



1 956 Provisional So uth American 
nates for the Schmidt telescope from Harrington ct 



Datu m. . ." . 
(1972]. 



al. 



aphi c sheet based on th e 1 
(d) |Walker Munoz I 
(e) iraf position in obsat 



the 1 924 International Reference El- 



1980[ l adopts the geodetic coordi- 
db.dat files: see e.g. |http:/ /tdc-| 

www.harvard.edu/iraf/ rvsao/bevcorr/ obsdb.html (/) MFC longitude was listed as 289° . 1941, with pcos(<ft ) — .86560 and 
psin(0) — -0.49980 (sec http://www.minorplanctccnter.net/iau/lists/ObsCodcs.html). (g) iHarring ton et al.| | |1972[ measured 
the astronomical position through timing the transits of stars with SAO catalog astrometry. They determined the geode- 
tic position through measuring and timing positions of the GOES-II satellite (which had a well-determined orbit, and strobe 
lights which flashed sequences of 7x 1.4 ms pulses every 4 s). The geodetic height is "above ellipsoid", (however which el- 
lipsoid is not explicitly mentioned (but given that the observations analyzed in 1971, it is likely WGS 66). (h) "O" is Tololo 
Control survey monument SE of Blanco 4-m. Values from AURA survey notes from Don Cassidy which triangulated the po- 
sition from Pcralillo and Pachon survey markers using theodolite observations, tied to elevations determined by Instituto Ge- 
ographic Militar (IGM). (i) NASA website (Space Geodesy Program): http://cddis.nasa.gov/sitc_cat/ccrr.html (j) LCOGT 
website: http:/ /lcogt.nct/sitc/ccrro-tololo (k) Gemini website: http: / /www. gcmini.edu/sciops/tclcscopcs-and-sitcs/locat ions 
(/) LSST website: http://www.lsst.org/lsst/science/summit_facilities (m) SOAR website: http://www.soartelescope.org/about- 
|soar /location-]"!] 



times, it may not be the position desired for other cal- 
culations (e.g. taking into account geocentric parallax, 
occultation calculations, etc.). 

3.3. 1973 Survey Elevations 

I list in Table [2] elevations for the concrete platforms 
for observatories on the Tololo plateau from a 1973 sur- 
vey map from the NOAO Engineering Department 20 . 
The concrete platforms for all of the observatories on 
the Tololo plateau have elevations ranging from 2209.60 
m (SMARTS/Yale 1.0-m) to 2210.50 m (Blanco 4-m), 

20 Contains note "Survey map drawn by E.W. Ross, checked & 
approved by Don Cassidy, Oct 12, 1971." 



i.e. less than AH = 1.1m among them. The elevations 
quoted in the survey are tied to topography from the 
"Fuerza Aerea de Chile - Servicio Aeoro Fotogrametrico 
(November 1964)" and "supplemented by plane table to- 
pography by AURA October 1966." The topography of 
the Tololo plateau was tied to elevations of neighboring 
peaks with coordinates provided by the IGM (including 
Cerro Pachon and Cerro Peralillo) via a primary survey 
monument called "Tololo", "Tololo Control ", or " O" on 
the old survey maps (discussed further in f 5.2.1). The 
Tololo survey monument was assigned elevation 2211.60 
m above sea level. Separate notes from this period list 
an elevation of 2724.60 m for an IGM survey monument 
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on Cerro Pachon, however it is not clear whether this 
monument corresponds to one of the modern day survey 
monuments on Pachon as the site has been developed 
considerably in the intervening four decades. 

3.4. 2008 Survey Elevations 

In Figure JT1 I show a scanned table from a 2008 survey 
of Cerro T ololo, Cerro Pachon, and the neighboring peak 
Cerro Morado (S of Tololo, NW of Pachon), completed 
by Juan Carlos Aravena Godoy for AURA. Cerro Morado 
has a Instituto Geographico Militar (IGM) survey mon- 
ument with well determined geodetic coordinates, which 
were provided to the surveyers by the Chilean Army 
(highlighted in red in Fig. n}. The coordinates are on 
the WGS 84 system (as adopted by SIRAS, the Sistema 
de Referenda Geocentrico para las Americas), and ele- 
vations are with respect to the GRS 80 ellipsoid. The 
difference in geocentric radii between the GRS 80 and 
WGS 84 ellipsoids at the position of the Blanco 4-m is 
26 /im (i.e. 10 -4,6 m). Hence, for all practical purposes, 
the differences between the GRS 80 ellipsoid adopted by 
SIRGAS (used in Chilean surveying) and the WGS 84 
ellipsoids, are completely negligible 21 . The 2008 survey 
also measured positions for 4 survey monuments near the 
Schmidt telescope on the Tololo plateau. The elevations 
ranged from 2241.175 to 2241.767 m, i.e. AH = 0.652 m, 
and H = 2241.430 m (elevations are with respect to the 
GRS 80 ellipsoid). This is probably representative of 
the GRS 80 heights of the concrete platforms for the ob- 
servatories on the Tololo plateau. The 1970s-era survey 
notes also list elevations for 4 pins distributed around the 
Schmidt telescope, with elevations ranging from 2209.333 
to 2209.985 m (AH = 0.592 m; H = 2209.650m), and 
given the description and close agreement in spread of 
elevations, these are likely to be the same pins as mea- 
sured in the 2008 survey. Both the Tololo Control and 
Schmidt survey monument measurements indicate that 
the 1973 survey maps can be converted to be geodetic 
heights with respect to the GRS 80 («WGS 84) ellip- 
soid by adding ^31 m. This difference is not far from 
the geoid undulations predicted for Tololo's position in 
recent EGM models (§5.2.2). These numbers seem to 
confirm that the 2008 survey elevations have not sub- 
tracted off geoid undulation, and they are simply with 
respect to the GRS 80 ellipsoid (hence, they are likely to 
have non-zero elevations at mean sea level). 

4. MEASUREMENTS 
GPS (Garmin): Positions were measured with a 
Garmin brand Dakota 20 GPS 22 , shown in Figure m This 
GPS was set to update at 1 second intervals, quoted po- 
sitions to 0".l precision and elevations to 1 m precision, 
and provided a regularly updated estimate of the posi- 
tion accuracy. Typical quoted accuracy in a 1 sec interval 
ranged from as good as 3 m (with unobstructed view and 

21 The only difference between the WGS 84 and GRS 80 ref- 
erence ellipsoid is due to the adopted flattening parameter (/) 
- in the 6th decimal place. So throughout, it is safe to assume 
that GRS 80 is equivalent to WGS 84. See also: NIMA Technical 
Report TR8350.2, "Department of Defense World Geodetic Sys- 
tem 1984, Its Definition and Relationships With Local Geodetic 
System s", 3rd edition, Amendment 1 (3 Jan 2000): http://earth- 
info.nga.mil/GandG/publications/tr8350.2/tr8350_2.htmlp 

2i Courtesy of Michael Warner (CTIO). 



SEC TOR CERRO PACHON 



Punto 


Latitud (s) 


Longitud (w) 


N o rt e (m) 


E s t e(m) 


Cota (m) 


Cerro Morado 


30°12 P 19,4380' 


70° 47' 46,691 3' 


6.657.090,474 


327.094,519 


2.175,292 


IGM Pachon 


30" 14' 26,7414" 


70° 44' 13,2744" 


6.653.259.605 


332.861,544 


2.748,353 


PE-4 


30° 14' 40,6827' 


70° 44' 57.5980" 


6.652.812.217 


331.683.277 


2.689.253 


WS-3 


30° 14' 32,0797" 


70° 44' 25,2738" 


6.653.090,346 


332.543,289 


2.750,504 


C° Tololo E-7 


30° 14' 42,81 93" 


70° 44' 47,9659" 


6.652.750.391 


331.941,765 


2.674,354 


SECTOR CERRO TOLOLO 










Punto 


Latitud (s) 


Longitud (w) 


N o rt e (m) 


E s tefm) 


Cota (m) 


Cerro Morado 


30" 12' 19,4380" 


70° 47' 46.691 3" 


6.657.090,474 


327.094,519 


2.175,292 


IGM C° Tololo 


30° 10' 11,7492" 


70° 48' 23,2709" 


6.661.006,347 


326.053,951 


2.242,493 


Tolo Este 


30" 10' 07,2465" 


70° 48' 12,781 1" 


6.661.149,421 


326.332,378 


2.183,515 


Pilole 1 


30" 10' 04,0150" 


70° 48' 18,3656" 


6.661.246,549 


326.181,406 


2.205,415 


SECTOR CERRO TOLOLO ESQ UENAS EDIFICIOS 






Punto 


Latitud (s) 


Longitud (w) 


N o rt e (m) 


E s t e(m) 


Cota (m) 


ED-1 


30° 10' 08,0850" 


70° 48' 21,5383" 


6.661.119,896 


326.098,513 


2.241,175 


ED-2 


30" 10' 07,5454" 


70° 48' 22,1604" 


6.661.136,244 


326.081,608 


2.241,229 


ED-3 


30" 10' 08,3004" 


70° 48' 22,8705' 


6.661.112,698 


326.062,978 


2.241,767 


ED-4 


30" 10' 08,6959" 


70° 48' 22,4386" 


6.661.100,706 


325.074,726 


2.241,550 



Fig. 1. — Geodetic positions for survey monuments on Cerro 
Tololo, Pachon, and Morado, from a report completed by Juan 
Carlos Aravena Godoy for AURA dated 28 October 2008. The 
coordinates are on the WGS 84 system (as adopted by SIRAS, the 
Sistema de Referencia Geocentrico para las Americas), and eleva- 
tions are with geodetic heights h respect to the GRS 80 ellipsoid 
(i.e. not "mean sea level", since no geoid has been subtracted). 
Measurements were taken via theodolite observations anchored to 
the Cerro Morado geodetic position (epoch 2002.0) provided by the 
Wilfredo Rubio Dias, Depto. de Datos y Calculo, Instituto Geo- 
graphico Militar, Ejerctio de Chile (dated 23 Oct 2008). Pins ED-1 
through ED-4 appear to correspond to positions near the Schmidt 
telescope on Tololo. 

several satellites acquired), and at times as bad as 20 m, 
close to the advertised 23 mean accuracy of ~15 m. 

Measurements were taken in different modes. As the 
GPS would not operate inside of domes due to obstructed 
view, I typically took two sequences of measurements at 
four corners of a given structure, and averaged the re- 
sults. For assymetric structures (e.g. Gemini, SOAR), 
I did some simple interpolation to estimate the central 
position. The GPS positions were very repeatable at the 
±0.1" level (only very rarely would subsequent measure- 
ments taken minutes apart differ by 0.2" ). Measurements 
were either taken (1) with the device sitting on the flat 
cement platforms outside of the domes, or the ground 
if there was no platform, or (2) with the device held 
in the hand, approximately 1 meter above the ground. 
None of the elevation measurements taken in this man- 
ner are quoted to better than 1 m precision. For GPS 
measurements of geodetic benchmarks and observatory 
sites currently lacking enclosures (e.g. KASI and T80- 
S), the GPS device was left sitting on the ground at the 
position, and coordinates were saved at 1 s or 5 s inter- 
vals for extended periods of time to the GPS's memory. 
As the GPS would often read out precisely the same co- 
ordinates for several seconds on end (then jump small 
amounts, presumably due to acquiring and losing indi- 
vidual satellites), it was decided to measure statistical 
moments only on unique coordinate sets in the time se- 
ries data. Especially long time series measurements were 
made of the benchmarks next to the SMARTS 1.5-m and 
NASA geodesy benchmark next to the SARA (former 
Lowell 24") telescope. 

GPS (iPhone): Auxiliary GPS measurements were 

23 http://www8.garmin.com/aboutGPS/ 
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TABLE 2 






Elevations for Observatory Platforms on 


1973 Survey Map 






Observatory 


H 


AH 




(m) 


(m) 



Blanco 4-m (158") 2210.50 -1.10 

SMARTS 1.5-m (60") 2210.99 -0.61 

SMARTS 1.0-m (40") 2209.60 -2.00 

SMARTS 0.9-m (36") 2210.49 -1.11 

Curtis Schmidt 0.6-m (24") 2210.00 -1.60 

No. 1 16" (later USNO) 2210.31 -1.29 

No. 2 16" (later CHASE) 2209.65 -1.95 

Tololo Control survey monument 2211.60 0.00 

Note. — H is elevation scale measured with re- 
spect to sea level, tied to topography established 
by Fucrza Acrca dc Chile - Servicio Aero Fotogra- 
mctrico, and supplemented by plane table topogra- 
phy measured by AURA. Later calculations show 
that the 1973 survey elevations appear to assume 
a gcoid within a meter or so of the EGM-84 model. 
AH is differential elevation with respect to Tololo 
Control survey monument. The mean elevation of 
the 7 platforms is 2210.22 m. Their mean height 
on Google Earth is 2202.7 m. Their mean height 
measured via GPS is 2217.6 m. 




Fig. 2. — The Garmin Dakota 20 GPS unit used in this study, 
sitting next to the benchmark adjacent to the southwest corner of 
the concrete platform for the SMARTS 1.5-m telescope. 

made with an iPhone 3s using Entel phone service, and 
using the application Compass. The iPhone measure- 
ments suffered two major limitations: sparse network 
coverage on the mountains, and the latitudes and lon- 
gitudes were listed at 1" precision. Elevation estimates 
were also made using the iPhone application Current El- 
evation 2 " 1 . As the quoted positions were at lower preci- 
sion than those provided by the Garmin GPS, and far 



|elevation/id429847648] 



http : // itunes.apple.com /us/app /current- 



fewer measurements were taken with the iPhone, I will 
not list these measurements. They do provide a consis- 
tency check on the Garmin GPS results, and indeed the 
iPhone coordinates appeared consistent with the Garmin 
results at the ~1" level, when available. 

Google Earth: I used Google Earth 5.1 downloaded 
from the Google website 25 and installed on a MacBook 
Pro running Mac OS X 10.6.8. Google Earth is a "virtual 
globe" that superposes satellite imagery over digital ele- 
vation model (DEM) data from the Shuttle Radar Topog- 
raphy Mission (SRTM). SRTM was a dual radar system 
that flew on Space Shuttle Endeavour in February 2000, 
producing a high resolution elevation model, and was 
a project lead by the Natio nal Geospatial-In telligence 
Agency (NGA) and NASA |Farr et al.||2007| . SRTM 
elevation data are on the WGS 84 coordinate system, 
adopt the EGM96 vertical datum (geoid), and for Chile 
the spatial resolution is 3" (~100 m). Hence, the Google 
Earth data from SRTM suffers from low resolution which 
may miss fine structure, and any leveling of sites over the 
past decade (e.g. LSST sites). Google Earth images of 
Cerro Tololo and Cerro Pachon are shown in Figures [3] 
and |4j respectively, with observatories and survey mon- 
uments marked. 

5. RESULTS 

Individual determinations of the geodetic positions for 
the observatories on Cerro Tololo and Cerro Pachon are 
tabulated in Tabled Additional measurements for other 
facilities on the two peaks are tabulated in Table [4] 

5.1. Comparison of Google Earth to GPS 

For approximately 30 structures (observatories, build- 
ings, and monuments) on Cerro Tololo and Cerro Pachon 
for which I measured both positions via GPS and 
Google Earth (GE), I find the following offsets in geode- 
tic latitude <fi and longitude A: 

<Pge ~ <Pgps = -0".10 ± 0".03(rms = 0".14) (2) 



25 http://www.google.com/earth/index.html 
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TABLE 3 
Measured Geodetic Positions 



Site 


Method 




4 


> 




A 


H 




















(m) 




Blanco 4-m 


GPS 


-30 


10 


10.73 


-70 


18 


23.52 


2213 


(fmr. 150" or 158") 


GE 


-30 


10 


10.73 


-70 


48 


23.32 


2204 


SMARTS 1.5-m 


GPS 


-30 


10 


09.53 


-70 


48 


24.49 


2218 


(former 60") 


GE 


-30 


10 


09.20 


-70 


48 


24.24 


2208 


SMARTS 1.0-m 


GPS 


-30 


10 


07.94 


-70 


48 


21.82 


2217 


(former Yale 40" ) 


GE 


-30 


10 


07.79 


-70 


48 


21.69 


2201 


SMARTS 0.9-m 


GPS 


-30 


10 


08.00 


-70 


48 


23.90 


2217 


(former 36") 


GE 


-30 


10 


07.69 


-70 


48 


23.67 


2203 


Curtis Schmidt 0.6-m 


GPS 


-30 


10 


08.59 


-70 


48 


22.62 


2219 


(former 24") 


GE 


-30 


10 


08.50 


-70 


48 


22.49 


2207 


Former UCAC 


GPS 


-30 


10 


06.93 


-70 


48 


22.85 


2222 


(former 16" #1) 


GE 


-30 


10 


06.87 


-70 


48 


22.65 


2200 


Former CHASE 


GPS 


-30 


10 


06.99 


-70 


48 


21.71 


2217 


(former 16" #2, MMWT) 


GE 


-30 


10 


06.84 


-70 


48 


21.44 


2196 


RASICAM 


GPS 


-30 


10 


07.99 


-70 


48 


25.10 


2214 




GE 


-30 


10 


07.96 


-70 


48 


24.97 


2200 


WHAM 


GPS 


-30 


10 


05.84 


-70 


48 


12.79 


2159 




GE 


-30 


10 


05.89 


-70 


48 


12.77 


2154 


KASI 


GPS 


-30 


10 


01.84 


-70 


48 


14.39 


2157 




GE 














2140 a 


LCOGT Stellan-A 


GPS 


-30 


10 


02.67 


-70 


48 


17.33 


2170 




GE 


-30 


10 


02.38 


-70 


48 


17.00 


2146 


LCOGT Stellan-B 


GPS 


-30 


10 


02.48 


-70 


48 


16.90 


2170 




GE 


-30 


10 


02.20 


-70 


48 


16.52 


2144 


LCOGT Stellan-C 


GPS 


-30 


10 


02.87 


-70 


48 


16.93 


2170 




GE 


-30 


10 


02.65 


-70 


48 


16.63 


2151 


SMARTS 1.3-m 


GPS 


-30 


10 


02.84 


-70 


48 


18.01 


2170 




GE 


-30 


10 


02.67 


-70 


48 


17.68 


2151 


PROMPT #1 


GE 


-30 


10 


03.42 


-70 


18 


18.77 


2162 


PROMPT #2 


GE 


-30 


10 


03.39 


-70 


18 


19.42 


2161 


PROMPT #3 


GE 


-30 


10 


03.07 


-70 


18 


18.71 


2157 


PROMPT #4 


GPS 


-30 


10 


03.59 


-70 


48 


19.40 


2178 




GE 


-30 


10 


03.42 


-70 


48 


19.09 


2162 


PROMPT #5 


GE 


-30 


10 


03.06 


-70 


18 


19.02 


2156 


PROMPT #6 


GE 


-30 


10 


03.72 


-70 


-18 


18.82 


2166 


PROMPT #7 


GPS 


-30 


10 


04.38 


-70 


48 


19.34 


2179 




GE 


-30 


10 


04.17 


-70 


48 


19.23 


2172 


GONG 


GPS 


-30 


10 


03.98 


-70 


48 


19.89 


2184 




GE 


-30 


10 


03.75 


-70 


48 


19.64 


2167 


SSI Airglow 


GPS 


-30 


10 


05.35 


-70 


48 


19.55 


2183 




GE 


-30 


10 


05.11 


-70 


48 


19.43 


2178 


SARA South 0.6-m 


GPS 


-30 


10 


19.73 


-70 


47 


57.11 


2123 


(Lowell 24" ) 


GE 


-30 


10 


19.61 


-70 


47 


57.00 


2112 


T80-South (site) 


GPS 


-30 


10 


04.31 


-70 


48 


20.48 


2187 


GE 














2175 a 


Gemini South 


GPS 


-30 


14 


26.54 


-70 


14 


12.13 


2722 




GE 


-30 


14 


26.70 


-70 


14 


11.84 


2711 


SOAR 


GPS 


-30 


14 


16.21 


-70 


14 


00.84 


2712 




GE 


-30 


14 


16.51 


-70 


14 


01.24 


2688 


LSST 8.4-m (site) 


GPS 


-30 


14 


40.68 


-70 


14 


57.90 


2652 


GE 














2633 a 


LSST Auxiliary 1.4-m (site) 


GPS 


-30 


14 


41.27 


-70 


14 


51.80 


2652 


GE 














2627 a 


Andes LIDAR Obs. (ALO) 


GPS 


-30 


15 


06.39 


-70 


14 


17.48 


2519 


GE 


-30 


15 


06.24 


-70 


14 


17.37 


2522 



Note. — GPS — Global Positioning System, measured using Garmin 
Dakota 20 unit. GE = Google Earth (26 March 2011 image of Tololo, 
11 Apr 2011 image of Pachon). Geodetic latitude <p, longitude A on 
WGS 84 system. Elevation H is orthomctric height with respect to gcoid. 
Google Earth adopts the EGM96 gcoid, however the adopted gcoid used 
by the Dakota 20 GPS is ambiguous. The KASI, T80-S, LSST, and LSST 
Support sites were undeveloped at the time of the latest Google Earth 
images (hence no coordinates could be determined visually), a — Google 
Earth elevation at the GPS position for this facility. The LSST and LSST 
Auxiliary sites were unleveled during the time (~2000) of the SRTM ele- 
vation mapping used by Google Earth, so arc not reliable. 
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Fig. 3. — Google Earth imagery of Cerro Tololo with observatories and survey monuments marked. 



TABLE 4 

Measured Geodetic Positions for Other Landmarks 



Site Method <f> A H 

o ' " o ' " (m) 



NASA Monument 7401 


GPS 


-30 10 20.919 


-70 48 00.096 2124.1 


(Site WZ, W of SaRa) 


GE 


-30 10 20.90 


-70 48 00.01 


2124 


Monument (SW of 1.5-m) GPS 


-30 10 09.76 


-70 48 24.93 


2214.0 




UK 


-30 10 09.67 


-70 48 24.78 


2206 


Tololo Control Monument 


GPS 
GE 


-30 10 11.75 


-70 48 23.25 


2222 


(SrD of Blanco) 






2201 


Round Office Building 
(Tololo) 


GPS 
GE 


-30 10 07.91 
-30 10 07.83 


-70 48 09.53 
-70 48 09.40 


2158 
2151 


AURA 20 Unit Dormitory GPS 


-30 15 04.38 
-30 15 04.34 


-70 44 20.80 
-70 44 20.75 


2502 
2506 


( Facnon ) 




New Dormitory 


GPS 


-30 15 07.54 


-70 44 23.64 


2503 


(PachonT 


GE 


-30 15 07.64 


-70 44 23.80 


2503 


New Dining Building 
( faction) 


GPS 
GE 


-30 15 07.26 
-30 15 07.22 


-70 44 24.01 
-70 44 23.97 


2503 
2500 


Old Dining Building 


GPS 


-30 14 29.25 


-70 44 17.76 


2699 


( faction) 


GE 


-30 14 29.28 


-70 44 17.90 


2692 


Communications Hut 
( Faction ) 


GPS 
GE 


-30 14 28.65 
-30 14 28.69 


-70 44 17.76 
-70 44 17.80 


2702 
2692 


CTIO Office Entrance 
(La Serena) 


GPS 
GE 


-29 55 01.27 
-29 55 01.37 


-71 14 31.01 
-71 14 30.97 


97 
90 


Gemini Office Entrance 
[La, Serena) 


GPS 
GE 


-29 54 59.50 
-29 54 59.65 


-71 14 30.31 
-71 14 30.44 


97 
94 


SOAR Office Entrance 
(La Serena ) 


GPS 
GE 


-29 55 02.35 
-29 55 02.38 


-71 14 31.85 
-71 14 31.21 


97 
88 



Note. — GPS — Global Positioning System, measured using Garmin 
Dakota 20 unit. GE — Google Earth. Geodetic latitude 0, longitude A on 
WGS 84 system. Elevation H is orthomctric height with respect to gcoid. 




LSST Support Telescope 



LSST 

?s tope *- 



AURA 20- 
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Googleearth 



Fig. 4. — Google Earth imagery of Cerro Pachon with observatories and buildings marked. 



X G e - A GPS = -0".14 ± 0".02(rms = 0".14) 
in units of length, this translates to: 



(3) 



i>GE - 4>gps = -3.1 ± 1.0 (rms = 4.3) m (4) 



Ag-E — ^gps = —3.8 ± 0.6 (rms = 3. 



(5) 



When comparing Google Earth imagery at different 
epochs, systematic shifts in latitude and longitude arc 
visible. Repeated measurements of the position for the 
NASA geodetic monument (just west of the SARA South 
observatory) on Google Earth at different dates (27 Feb 
2006, 3 Apr 2010, 26 Mar 2011) show systematic epoch- 
to-epoch shifts at the ±0". 08-0". 13 level. Hence, for the 
purposes of deriving "best" geodetic positions, I correct 
the Google Earth positions to take into account their 
systematic difference with respect to GPS positions. 

5.2. NASA Geodetic Monument 

The NASA Crustal Dynamics Project established a 
cluster of observing monuments near the SARA South 
Observatory (former Lowell 24") for satellite laser rang- 
ing measureme nts (the technique is discussed in e.g. Tap- 
ley et al.|198 5). The monuments are collectively referred 



to as "site number 892" 26 , but the primary monument is 
labeled number 7401, and it occasionally appears by this 
number in the geodesy literature. The monument disk 
is about a meter north of the center of a 25 foot square 
pad easily visible on Google Earth (and easily visible just 
south of the SARA telescope access road), and shown in 
Figure 5} The site hosted three campaigns using Trans- 
portable Laser Ranging Systems (T LRS) between 1984 
and 1991 (the program is discussed in Allenby|1984[ ), tak- 
ing range measurements to the Laser Geodynamics Satel- 
lite (LAGEOS). Precise geodetic coordinates for monu- 
ment 7401 are provided in the SGAPO online archives as 
[ latitude south 30° 10' 20.86692" , longitude west 70° 48'| 
|00.15364"[ elevation 2123.090 m, height above ellipsoid 
2155.748 in (dated 23 April 1990). The quoted ellipsoid 
assumed an equatorial radius of 6378137 m and flatten- 
ing f = 1/298.255, i.e. an identical radius, but slightly 
different flattening compared to WGS 84. 

The monument provides a useful check on the types 
of "elevations" that are being reported. Both GPS and 
Google Earth measured an elevation of 2124 m for the 
NASA monument, i.e. only 1 m above the NASA geode- 
tic value. This measurement alone is strongly suggestive 
that both the GPS and Google Earth measurements take 
into account a geoid model, and are not simply measured 
from the WGS 84 reference ellipsoid. Unfortunately, the 



1>G 



Data 



provided 



at the website 
http://cddis.nasa.gov/site_cat/cerr.html maintained by the 
Space Geodesy and Altimetry .Projects Office (SGAPO), edited 
by Mark Bryant and Carey Noll, dated March 1993. 
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Fig. 5. — Platform for monument 7401 (site #892) situated just 
west of the SARA South Observatory (former Lowell 24"). The 
benchmark was built for the NASA Crustal Dynamics Project in 
the early 1980's, and satellite laser ranging measurements were 
taken here between 1984 and 1991. The actual benchmark is 
marked by a metal disk within the central dark square on the plat- 
form. Several auxiliary benchmarks labeled with 7401 and letters 
are visible scattered across the platform. The SARA South Obser- 
vatory lies to the right of the image down the service road that the 
car is parked on. 



agreement between the GPS and Google Earth elevations 
for the monument does not provide a useful explanation 
for the ^11 m systematic offset between the two as in- 
ferred from averages of elevation measures for ~30 other 
sites. The source of the discrepancy at other sites is un- 
clear. Is standing next to buildings biasing the GPS ele- 
vations? Is small scale elevation structure not taken into 
account by the smoothed SRTM elevation data affecting 
the Google Earth elevations? 

Preliminary es timates of the motio n of monument 7401 
were reported in Smith et al. ( 1994 ) , however the errors 
were large. The measured motion during 1984-1991 was 
35.9 mm/yr towards azimuth angle 30°, however the er- 
ror ellipsoid was 24.0 mm/yr x 10.7 mm/yr with the 
long axis oriented towards azimuth angle -5°. Hence, the 
measured motion was statistically negligible. The NASA 
geodetic coordinates are valuable, however, as they pro- 
vide not only a well-calibrated position for comparison 
with the GPS and Google Earth positions, but a poten- 
tial first epoch for estimating the motion of Cerro Tololo 
with respect to the WGS 84 terrestrial coordinate sys- 
tem. 



5.2.1. Tololo Control Monument 



A survey monument called "Tololo" , "Tololo Control" , 
|or "Q"|o n the 1960s/1970s-era notes and survey maps of 
Gerro Tololo lies just southeast of the Blanco 4-m dome 
structure, just outside of the guard rail girding the gravel 
(see Fig. |6| . In Table [I>1 I list several estimates of the 
elevation of the Tololo "Control survey monument] An 
elevation of 2211.60 m above sea level is listed on a sur- 
veying note from the Kitt Peak National Observatory 
Engineering Dept. This benchmark was used in conjunc- 
tion with monuments on other nearby peaks (Pachon, 
Morado, etc.) for tying in the position of Tololo with 
the cartographic grid used by the Instituto Geographico 
Militar de Chile and Chilean Air Force in the 1960s (the 




FlG. 6. — "Tololo Control" survey monument on the Cerro Tololo 
peak, just southeast of the Blanco 4-m telescope. This was the pri- 
mary benchmark for Tololo in the surveys of the mid-1960s, and its 
elevation is within ~l-2 m of the platforms for all of the telescopes 
on the Tololo plateau. The 2008 survey by Juan Carlos Aravena 
Godoy estimated a geodetic position of </> = -30°10'11.7492" , A = 
-70°48'23.2709" at elevation 2242.493 m above the GRS 80 ellip- 
soid (i.e. not taking into account a geoid model, not with respect 
to mean sea level). 

Chile Plane Coordinate System) 27 . The elevation of this 
monument defined the elevation scale for other structures 
built on the Tololo plateau, as marked in the 1973 survey 
map (see Table [2~|. GPS measurement of the elevation 
of the |Tololo m onument yielded 2222 m. The 2008 sur- 
vey discussed in §3.4| measured a geodetic elevation of 
2242.493 m with respect to the GRS 80 ellipsoid (which 
is within tens of micrometers of the WGS 84 ellipsoid). 
The monument is not obvious on Google Earth imagery, 
however the elevation in the vicinity of the monument is 
2201 m. 



The 2008 survey measurement for the Tololo Con- 
trol survey marker is the most accurate WGS 84 co- 
ordinate position, and fortunately its datum is unam- 
biguous (GRS 80 « WGS 84). The GPS elevation is 
20 m lower than the 2008 geodetic survey elevation. The 
Google Earth elevation is ^10 m lower yet, however 
Google Earth measures elevations with respect to the 
EGM-96 geoid, which accounts for most of the ^41 m 
discrepancy. The level of disagreement among the eleva- 
tions is somewhat surprising, however, given the good 
agreement for the elevations measured for the NASA 
monument near the SARA South observatory. 

5.2.2. Treatment of Elevation Data 

Multiple series of GPS measurements were taken of 
a few sites on different days. From comparing average 
GPS elevation measurements taken on different days, it 
is clear that significant systematic errors in the GPS el- 
evations are present at the ±5 m rms level (presumably 
from having data from different satellites participate in 
the solution). This precludes using the standard error of 
the mean for measurements taken during a short interval 
on a given day as a useful measure of the uncertainties. 

Inter-comparison of the GPS and Google Earth eleva- 
tions reveals a systematic offset: 

27 The monument disk states "Instituto Geographico Militar de 
Chile - destruccion penada por la ley" . 



Coordinates for Cerro Tololo and Cerro Pachon 
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TABLE 5 

Elevations for Tololo Control Monument 



Source Elevation 
(m) 




Notes 




1973 survey 2211.60 m 
2008 survey 2242.493 m 

Google Earth 2201 m 
GPS 2222 m 


"mean sea level" 

w/r GRS 80 ellipsoid (no geoid) 

w/r WGS 84 ellipsoid + EGM96 geoid 

w/r WGS 84 ellipsoid + unknown geoid 



< H Google - H GPS >= -11.5 ± 2.0 



(6) 

The rms scatter was ±8.6 meters. The source of the 
discrepancy is unclear. Two likely possibilities are that 
either (1) the 3" resolution of the Google Earth elevation 
data from SRTM has smoothed over the high elevation 
points (where many of the telescopes are located), or (2) 
there may be a difference in the geoid adopted in the GPS 
calculations and that used by Google Earth (EGM96). 
Both the GPS and Google Earth elevations are clearly 
limited in their accuracy. The GPS measurements often 
have epoch-to-epoch jumps and the datum is ambiguous 
(at the time of writing, there is conflicting literature on 
the web regarding this; whether Garmin GPSs measure 
with respect to the ellipsoid, or whether they include a 
geoid model). The Google Earth maps have excellent 
imagery resolution, but the elevation maps are limited 
to the SRTM resolution (~3"). Unfortunately, elevations 
taken both with the handheld GPS unit and via Google 
Earth must both be taken with a grain of salt. The best 
elevations that we have available are those from the 2008 
survey and those for a NASA geodetic monument. 

For useful comparison of the elevations, it was de- 
cided to place all measurements on the same elevation 
system defined by the 2008 theodolite survey (which in- 
cludes very accurate elevations measured with respect 
to the GRS 80 [«WGS 84] ellipsoid). Through inter- 
comparison of the measured elevations, I decided to cor- 
rect GPS, Google Earth, and 1973 survey elevations to 
the geodetic elevation system defined by the 2008 survey 
(close enough to WGS 84 to be practically considered 
WGS 84): 



h\VGS 



Hi 



Hgps + (28 ± 3) m 



(7) 



h\VGS 84 — H Google Earth 

+ N(EGM96) (8) 



h\YGS 84 — #1973 Survey + (31.3 ± 0.4) m (9) 

The difference between the Google Earth elevations 
and the 2008 survey geodetic heights h are close enough 
to the EGM-96 geoid undulations N (38 ± 3 m; see ac- 
tual EGM-96 N values in Table 6), that they appear 
to provide useful elevations H . The difference between 
the GPS elevations and 2008 survey heights do not ap- 
pear to correspond to either zero (if the GPS elevations 
were actually geodetic heights), or any of the widely used 
geoid undulations. Hence, I simply add a constant to the 
GPS elevations and place them on the WGS 84 geode- 
tic height scale. The difference between the 2008 survey 
geodetic heights and the 1973 survey elevations (31.3 m) 
is remarkably similar to the geoid undulation heights (§2; 



especially for the WGS 84 ellipsoid) , suggesting that the 
zero point for the 1973 survey elevations is indeed within 
a few meters of recent geoid models. 

6. CONCLUSIONS 

Best estimates for the geodetic and geocentric positions 
of the observatories on Cerro Tololo and Cerro Pachon 
are listed in Table [6j Final elevations taking into ac- 
count the offsets in elevations discussed in the last section 
are also included, along with geoid undulation heights N 
based on the EGM-96 model. 

It appears that using Google Earth one can derive 
geodetic WGS 84 coordinates for observatories and struc- 
tures to approximately ±0".15 (5 m) accuracy in lati- 
tude and longitude without resorting to more accurate 
techniques (i.e. GPS, surveying). Google Earth shows 
systematic shifts in position between imagery epochs at 
the ~±0".l level (~3 m). Positions measured on the 
most recent Google Earth imagery of Cerro Tololo and 
Pachon reveal slight systematic offsets with respect to 
GPS-derived positions at the ±0".l (3 m) level, which 
can be corrected for if necessary. The rms scatter in 
positions derived via GPS and Google Earth is approxi- 
mately 0".14 (4 m). Treating the GPS and Google Earth 
positions as independent estimates, and correcting the 
Google Earth positions for small systematic differences at 
the 0".l (3 m) level, one can derive final mean positions 
for structures on Cerro Tololo and Pachon with abso- 
lute accuracy ±0".l (3 m). One can obviously do better 
if needed through long integrations of GPS determina- 
tions, or theodolite observations tied to monuments with 
well-constrained geodetic positions, however this would 
obviously require more time and money (which always 
seem to be in short supply). 

While astronomical coordinates have not been deter- 
mined for all of the observatories on Cerro Tololo and 
Cerro Pachon (an arduous task), one can derive approx- 
imate astronomical coordinates by adding an offset de- 
rived from Harrington et al.'s (1972) observations: 



" astronomical 



"geodetic 



X. 



astronomical 



■><!< I /■ 



12". 7 ± 1".5 



-31" .6 ± 2".0 



(10) 



(11) 



However this is only a rough approximation as no 
doubt the vertical deflection due to the gravity field 
varies subtly over the two peaks. 

Besides the obvious (and explainable) disagreement be- 
tween th e ge odetic and astronomical positions for Cerro 
Tololo ({3.2 1, there was at least one other surprise in 
comparing the final positions in Table [6] with the previ- 
ously published positions in Table [l] ^n often quoted 
position for LSST listed in two recent LSST documents 



12 



Mamajck 



(Ivezic et al. 2008a|b ) appears to be in error by 9.4 km 
[Acj) = 260" .6, AA = 60.0"), and by ~500 m in elevation 
(2123 m listed vs. ^2600 m mea sured). The position 



listed in |Ivezic et aL] ( |2008a[ ) and |Ivezic et al] ( |2008b| ) 
appears to correspond most closely to the IN AS A survey 
monument on Cerro Tololo (rather than Pachon). The 
position for the LCOGT observatories on Tololo on the 
LCOGT website appear to be in error by ^0.2 km. Our 
position for Gemini South agrees with the value on their 
website 28 to better than <0".04 for each axis (within ~1 
m), and our GPS estimate of the elevation (2722 m) is 
identical to their value. 

After finding that some of the published positions for 
observatories on Cerro Tololo and Cerro Pachon were er- 
roneous, it may be worth the time to double-check the 
published positions of other observatories (especially if 
those observatories pre-date the GPS era). Using erro- 



sitcs/locations 



|http://www.gemini.edu/sciops/telescopes-and-| 



neous positions for observatories can lead to unneces- 
sary systematic errors in various position-dependent as- 
tronomical calculations, which are inexcusable in the era 
of GPS and Google Earth. 



The author thanks Jason Wright for pointing out the 
discrepancy between the published position for CTIO 
and that determined with Google Earth, which triggered 
this investigation. I thank Michael Warner for lending 
the author his Garmin GPS. The author also thanks En- 
rique Figueroa for sharing the 2008 survey report from 
Juan Carlos Aravena Godoy, and thanks to Malcolm 
Smith, Oscar Saa, Tim Abbott, Alistair Walker, Chris 
Smith, Nicole van der Bliek, Kadur Flores, Jeff Barr, 
and Dan Phillips, for conversations on the history and 
mythology of Cerro Tololo and Cerro Pachon. 
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Changes 

4 Oct 2012: 
9 Oct 2012: 
in Tables 1, 
changed. 



http://xxx.lanl.gov/abs/1210.1616j 

ATided Google Maps hyperlinks for all coordinates listed 



Posted to arXiv 

Minor edits. Added hyperlinks for all websites 
4, and 6. Table 6: Added footnote on LSST floor elevation, and added tablenotemarks 



No numbers were 



Coordinates for Cerro Tololo and Cerro Pachon 



TABLE 6 

Geodetic and Geocentric Coordinates for Observatories on 
Cerro Tololo and Cerro Pachon 



Site 
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ietic 


A 


geodetic 


Ref. 


h 


N 


H 
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O 








o 










(m) 
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(m) 
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-70 


48 


19 


56 


1 
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04 
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48 
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2 
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34.6 
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58 


SARA South O.e 


-m 


-30 


10 


19 


72 
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47 


57 


13 


1 


2151 c ' d 


34.7 
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-30 


00 


18 


93 



(Jerro Hachon 



Gemini South 8.2-m 
SOAR 4.i-m 



LSST 8.4-m (site) 



LSST Aux. 1.4-m (site) 



Andes LIDAR Obs. (ALO) 



-30 14 26.67 
-30 14 16.41 
-30 14 40.68 
-30 14 41.27 
-30 15 06.37 



70 44 12.06 1 
70 44 01.11 1 
70 44 57.90 2 
70 44 51.80 2 
70 44 17.50 1 



2750 c 
2748 e 
2647 / 
2647 / 
2552 c >< 



35.1 2715 c 
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-30 04 25.12 
-30 04 14.89 
-30 04 39.07 
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Note. — Coordinate reference codes: (1) average of GPS and corrected Google Earth positions, (2) 
GPS position, (3) corrected Google Earth position. Geodetic latitude 0, longitude A on WGS 84 system. 
Geodetic elevation h is with respect to the WGS 84 ellipsoid. Elevation H is orthomctric height with 
respect to gcoid (i.e. above "mean sea level"; H — h - N). N is gcoid undulation according to EGM-96 

model, i.e. it is the height of t he geoid above the WGS 84 reference ellipsoid. Geocentric latitude (ft 
was calculated following §4 of |Seidelmann| l |1992| . Total uncertainties in latitude and longitude are 
approximately ±0".10. Elevations arc ground level at the site, not the elevation of the telescope. 
a Ground elevations of facilities on Tololo plateau arc likely accurate to better than ±1 m absolute 
(calculated using 2008 survey geodetic height for Tololo Control survey monument and differential 
heights from 1973 survey). 

According to building plans for the 1.5-m, the center of its dome is 10.3 m above ground level, so the 
1.5-m telescope is at h ^ 2252 m and elevation H ~ 2218 m. 

c Elevations of other Tololo facilities, Gemini, and ALO arc a blend of corrected GPS and Google 
Earth elevations, but arc likely accurate to better than ±5 m absolute. 

^ The elevations listed for facilities on the slopes of Cerro Tololo are commensurate with that of 
a survey monument labeled "Pilote 1" (/iG-RS 80 — 2205.415 m ~ hyvGS 84) situated between the 
SMARTS 1.3-m and PROMPT cluster. 

c For the elevation of SOAR, I adopt the 2008 survey data for Pachon monument next to SOAR 
telescope. 

^ Elevations for LSST and LSST auxiliary telescope are adopted from LSST website 
http://www.lsst.org/lsst/scicncc/summit_facilitics which were based on cartography tied to the 2008 
survey. According to current LSST plans, the pier floor under the telescope will be at h = 2662.75 m. 



